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^ ■ ABSTRACT 

j- \ We present the morphological analysis based on HST-NIC2 (0.075 arcsec/pixel) images 

(SJ ■ in the F160W filter of a sample of 9 massive field (> 10^^ Mq) galaxies spectroscopi- 

' cally classified as early-types at 1.2 < z < 1.7. Our analysis shows that all of them are 

bulge dominated systems. In particular, 6 of them are well fitted by a de Vaucouleurs 
i profile {n = 4) suggesting that they can be considered pure elliptical galaxies. The 

' remaining 3 galaxies are better fitted by a Sersic profile with index 1.9 < n/it < 2.3 

r~| . suggesting that a disk-like component could contribute up to 30% to the total light 

Q^' of these galaxies. We derived the effective radius Re and the mean surface brightness 

(/x)e within Re of our galaxies and we compared them with those of early- types at 
lower redshifts. We find that the surface brightness (^)e of our galaxies should get 
^ fainter by 2.5 mag from z ~ 1.5 to z ~ to match the surface brightness of the local 

^ \ ellipticals with comparable Re, i.e. the local Kormendy relation. Luminosity evolution 

without morphological changes can only explain half of this effect, as the maximum 
dimming expected for an elliptical galaxy is '^1.6 mag in this redshift range. Thus, 
[v> , other parameters, possibly structural, may undergo evolution and play an important 

j_j ■ role in reconciling models and observations. Hypothesizing an evolution of the effec- 

' tive radius of galaxies we find that Re should increase by a factor 1.5 from z ^ 1.5 to 

z ~ 0. 

Key words: Galaxies: evolution; Galaxies: elliptical and lenticular, cD; Galaxies: 
formation. 
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1 INTRODUCTION 

The working out of the hierarchical clustering paradigm 
paved the way to the development of various models of 
galaxy formation. All these models are based on the assump- 
tion that gravity is the main driver shaping the structures, 
both on galactic and larger scales. Observational evidences 
sustaining this scenario, in which structures form by subse- 
quent mergers of dark matter haloes, come from the WMAP 
data which constrain the fundamental parameters of mod- 
ern cosmology (f2o, the Hubble constant and the amount 
of dark energy) to better then a few per cent (e.g. Spergel 
et al. 2003). However, while the general build-up of cos- 
mic structures seems to be well described, the assembly of 



* E-mail: paolo.saracco@brera.inaf.it 



the baryonic mass on galactic scales still represents a weak 
point and it does still not find a satisfactory description in 
these galaxy formation models. In particular, the assembly 
of high- mass spheroids is a key issue of the current models of 
galaxy formation since they still miss the process accounting 
for their growth. 

As recently pointed out by Nagamine et al. (2005a) , one 
of the difficulties of hierarchical models consists indeed in ac- 
counting for the population of massive (M,t,, > 10" Mo) 
early-type galaxies seen at high z. Indeed, while they show 
that high-mass optically selected (blue) galaxies at z ~ 3, 
i.e. the Ly-break galaxies, can be easily recovered in the 
framework of hierarchical models (Nagamine et al. 2QQ5b), 
the number of massive early-type (red and dead) galaxies 
seen at z > 1 hardly fit the predictions of these models (see 
also Somerville et al. 2004). A recent important result has 
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been, indeed, the discovery and the spectroscopic identifica- 
tion of a significant population of massive early-type galax- 
ies at 1.2 < z < 1.7 (McCarthy et al. 2004; Glazebrook 
et al. 2004; Saracco et al. 2005). Their number density is 
comparable to the number density of massive early-types 
seen in the local Universe showing that no strong evolution 
has taken place up to x ~ 1.5 (Saracco et al. 2005). On 
the other hand, their spectrophotometric properties reveal 
a spread in the mean age of their stellar population sug- 
gesting a diffcreut star formation history and/or a different 
redshift of formation (Longhetti et al. 2005). The discovery 
of early-type galaxies at 2; ~ 2 (Cimatti et al. 2004; Daddi 
et al. 2005) shows that early-types exist at least up to ~ 2. 

Another important observational evidence is the mass- 
dependent evolution which the galaxies undergo apparently 
(Gavazzi 1993; Gavazzi et al. 1996), the so called "down- 
sizing" (Cowie at al. 1996). This observed behaviour con- 
strains the build up of the stellar mass in high-mass early- 
type galaxies earlier and over a shorter interval than in the 
lower mass ones. This behaviour has been recently confirmed 
observationally both at low and at high redshift through a 
variety of analysis: from the abundance ratios of local early- 
type galaxies (Thomas ct al. 2005), from the evolution of 
the fundamental plane (Treu et al.2005), of the specific star 
formation rate (SSFR; Feulner et al. 2005) and of the stel- 
lar mass function (Drory et al. 2005; Bundy et al. 2005) of 
galaxies, from the clustering of distant red galaxies (Grazian 
et al. 2006), from the optical LF of galaxies at 2 < 1 (Cimatti 
et al. 2006) and from the evolution of the near-IR LF of 
galaxies to 2: ~ 3 (Saracco et al. 2006; Caputi et al. 2006). 

These pieces of evidence arc apparently at variance with 
the expectations of hierarchical models which predict a z > 2 
Universe mostly populated by star-forming disks, irregular 
galaxies and merging systems and in which the local massive 
spheroids complete their assembly last, possibly at 2; < 1. 

To further complicate this picture, evidence for a pos- 
sible lack of big early-type galaxies at 2: ~ 1 and beyond are 
coming out (e.g. Fasano et al. 1998; Waddington ct al. 2002; 
Daddi et al. 2005; Trujillo et al. 2006a, 2006b; Cassata et al. 
2005). At first glance, this could be interpreted as a prove of 
the merging process responsible of the growth of local high- 
mass ellipticals. On the other hand, the apparently small 
early-types seen at high-redshift seem to be characterized by 
an effective surface brightness brighter than the one of local 
counterparts even considering luminosity evolution. More- 
over, their stellar masses exceed those of local galaxies with 
comparable size. Thus, they seem to be more compact then 
their local counterparts. It should be noted, however, that 
most of these results arc based on HST optical observations 
sampling the blue and UV rest-frame emission of the galax- 
ies, particularly sensitive to morphological k-correction and 
star formation episodes, and/or on seeing limited ground- 
based observations. These features could affect the estimate 
of the effective radius of high-z early-type galaxies. 

In this paper, we present the morphological analy- 
sis based on high-resolution (FWHM^O.l arcsec) HST- 
NICMOS observations in the F160W filter (A ~ 1.6 urn) 
of a sample of 9 high-mass early-type galaxies with spec- 
troscopic confirmation at 1.2 < z < 1.7. The filter used 
samples the rest-frame R-band at the redshift of our galax- 
ies. The aim of this work is to establish whether the surface 
brightness vs. effective radius relation of early-type galaxies 



Table 1. Properties of the massive early- type galaxies of our 
sample. Details on all the reported parameters can be found in 
Longhetti et al. (2005). 
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at 2 > 1 can be accounted for by their luminosity evolution 
or whether they are in fact more compact than local early- 
types without the uncertainties related to the optical and 
ground based observations. 

The paper is organized as follows: in §2 we describe the 
sample and the HST observations. In §3 we describe the 
profile fitting procedure used, we derive the morphological 
parameters of our sample of early-types and we describe the 
simulations performed to assess the robustness of our fitting. 
In §4 we report the results of our morphological analysis. In 
§5 we discuss the Kormendy relation, we compare our results 
with those obtained at lower redshift by other authors, we 
constrain the evolution that early-types must undergo at 
2 < 1.5 and we propose a model accounting for the observed 
evolution. Summary and conclusions are given in §6. 

Throughout this paper we assume Vega magnitudes, 
Ho=70 km s"^ Mpc"^ Ho = 0.3 and Ao = 0.7 flat cos- 
mology. 



2 SAMPLE AND HST-NICMOS 
OBSERVATIONS 

The sample is composed of 9 out of the 10 field galaxies 
spectroscopically classified as early-types at 1.2 < 2 < 1.7 
{{z) = 1.4, Longhetti et al. 2005). They result from a near- 
IR spectroscopic follow-up of a complete sample of bright 
(K< 18.5) EROs (R-K> 5.3) selected from the Munich 
Near-IR Cluster Survey (MUNICS, Drory et al. 2001) that 
provides optical (B, V, R, I) and near IR (J and K') pho- 
tometry. The low resolution spectroscopic observations were 
carried out with the prism dispersor AMICI at the near- 
IR camera NICS of the TNG (Tcloscopio Nazionale Galileo, 
La Palma, Canary Islands). The spectroscopic and the pho- 
tometric data show that they have stellar masses greater 
than 10^^ Mq and that they account for more than 70% of 
the number and the stellar mass density of local early-type 
galaxies with comparable masses (Saracco ct al. 2005). The 
stellar population is 3-5 Gyr old in 6 of them and 1-2 Gyr 
old in the remaining 4 galaxies implying minimum formation 
redshifts 2/ > 4 and 2/ > 2 respectively (Longhetti et al. 
2005). The properties of the 9 early- type galaxies (absolute 
K-band magnitude, stellar mass and age) are summarized 
in Tab.l. 

The near-IR HST observations were carried out dur- 
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ing Cycle 14 with the NICMOS-NIC2 camera (0.075 arc- 
sec/pixel) in the F160W filter (A ~ 16030 A). Consider- 
ing the redshift of our galaxies, the selected filter samples 
their R-band rest-frame continuum. Each galaxy has been 
observed for 2 orbits (~ 4800 seconds exposure) with the 
exception of the brightest one (S2F5_109) for which we al- 
located a single orbit (2400 sec). The observations consist 
of 8 dithered exposures of about 300 seconds each for each 
orbit. The MULTIACCUM Read Out Mode STEP16 with 
NSAMP=25 and a dither size of 4 arcsec with a spiral dither 
pattern have been adopted. 

The HST calibrated images used for our analysis are 
those provided by the HST data archive resulting from the 
standard automatic reduction procedures CALNICA and 
CALNICB. The limiting surface brightness reached in our 
images is hfigo — 27.77 mag/pix (29.11 AB mag/pix) 
corresponding to /Ufiso — 22.14 mag/arcsec^ (23.48 AB 
mag/arcsec'^). This allows us to sample the light profile of 
our galaxies at a S/N>3 out to more than 3 effective radii. 

3 MORPHOLOGICAL PARAMETERS 

The measure of the morphological parameters of the galax- 
ies, such as their effective radius (re in arcsec or Re in kpc) 
and the mean effective surface brightness within Re {{fije in 
mag/arcsec^), has been performed by means of the Galf it 
software (v. 2.0.3; Peng et al. 2002). Galf it builds a bi- 
dimensional model of the image of the galaxies on the basis 
of an analytic function that can be chosen by the user. It 
takes into account the position angle (PA) of the major axis 
a of the profile and the ratio a/b between major and mi- 
nor axis. The bi-dimensional model is convolved with the 
Point Spread Function (PSF) of the observations before be- 
ing compared with the observed image. Minimization of the 
differences between models and observations gives the best 
fitting parameters which have been left varying in the pro- 
cedure. 

We adopted a single component to fit the light distri- 
bution of the galaxies in our sample. We performed the fits 
adopting in separate steps the de Vaucouleurs (1948), the 
exponential and the Sersic (1968) profiles. The analytic ex- 
pression of the adopted profiles is 

7(r) = hexp{-b4{r/r,)^/" - 1]} (1) 

where /(r) is the surface brightness measured at the dis- 
tance r from the centre of the galaxy in flux units/arcsec"^, 
Je is the surface brightness measured at r = re where re 
[arcsec] and b„ is a normalisation factor depending on the 
exponent n chosen for the fit. The exponent n is a free pa- 
rameter for the Sersic profile, while it is n = 1 and n = 4 
for the exponential and de Vaucouleurs profiles respectively. 
The other free parameters in the adopted fits arc the total 
magnitude (magcoj), the semi-major axis Oe of the projected 
elliptical isophote containing half of the total light, the ax- 
ial ratio b/a and the position angle of the major axis. The 
circularised effectiv e ra dius is related to Oe through the 
equation re = a^^/bja. The mean effective surface bright- 
ness (/Lt)e can be computed considering half of the total flux 
of the galaxy and the area included within the effective ra- 
dius 



The magnitude in the F160W filter used to derive (/i)e is 
the total magnitude magoai derived by Galf it. 

A different PSF has been generated for each galaxy by 
means of the Tiny Tim^ (v. 6.3) software package (Krist 
1995; Krist J. & Hook R. 2004). Since the final image of each 
galaxy is the result of the co-adding of 8 shifted images, the 
relevant PSF has been constructed by averaging the 8 PSFs 
generated at each target position. No sub-sampling has been 
adopted and the best fitting Spectral Energy Distribution 
(SED) of each galaxy (Saracco et al. 2003, 2005; Longliotti 
et al. 2005) has been used to simulate the spectral shape 
of the light entering the F160W filter. We verified that the 
results of the fitting do not depend on the specific choice of 
the Tiny Tim parameters mentioned above. In particular, we 
find that the best fitting parameters vary less than 2% (e.g. 
A(re) < 0.01") for different choices of SED, sub-sampling 
factors and PSF position. 



3.1 Profile fitting 

We first assess whether the light distribution of our galax- 
ies is bulge-dominated, as expected from our spectroscopic 
classification, or disk-dominated. To this end we compared 
the goodness of fit obtained with the de Vaucouleurs (n = 4, 
bulge) profile and the exponential (n = 1, disk) one making 
use of the relevant x^- The fitting has been performed over 
a box of 6x6 arcsec centered on the centroid of each galaxy 
after having checked for the stability of the results with re- 
spect to the size of the box. As initial guesses for the cen- 
troid, for the PA and for the total magnitude of the galaxies 
we used the SExtractor output. We did not keep the SEx- 
tractor MAG_BEST value as the actual total magnitude of 
the galaxies since it turned out to be less accurate than the 
one resulting from Galf it, as we verified through simula- 
tions (see next sub-section). Neighbouring sources and bad 
columns resulting from the mosaicking have been masked 
out from the model fitting. The root moan square imago pro- 
duced by the standard HST-NICMOS reduction procedure 
has been used to give the relative weights to the sky pixels 
in the fitting. We find that the de Vaucouleurs profile always 
provides a better fitting than the exponential one for all the 
galaxies. In particular, the n = 1 profile provides a fitting 
significantly worse (A(x^) > 0.3, d.o.f.> 100) than the de 
Vaucouleurs one for 8 early-type galaxies while, for galaxy 
S2F1_443, the difference is not significant (A(x^) ~ 0.01). 
It is interesting to note that S2F1_443 is the highest red- 
shift galaxy of our spectroscopic sample of early-types. Thus, 
both the surface brightness dimming and the angular scale 
are larger than for the others. Besides this, this galaxy is 
"peculiar" with respect to the other galaxies since it is X- 
ray emitting, probably hosting a type 2 QSO (Severgnini et 
al. 2005). 

From this first comparison, we cam exclude that a disk 
component is the dominant feature of the light distribution 
of our early-type galaxies, in agreement with their spec- 
trophotometric properties. In Fig. 1 the F160W-band im- 
age, the best-fitting de Vaucouleurs model and the residual 
of the fitting are shown for each galaxy. 



{Ai)e = F\mW{mag) -|- Uog{re) + 2.5Zog(27r) 



(2) 
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Figure 1. Each column shows the Galf it input and output for the 9 early-type galaxies: NIC2 F160W-band image of the galaxy (upper 
panel), best- fitting de Vaucouleurs model profile (middle panel) and residual image (lower panel) obtained by subtracting the model 
from the image. Galaxies are (from top-left to bottom-right): S2F5.109, S7F5.254, S2F1_357, S2F1_389, S2F1_511, S2F1_142, S7F5.45, 
S2F1_633, S2F1_443. Each image is 6x6 arcsec. 



We then compared the Sersic profile to the de Vau- 
couleurs profile to verify whether the addition of the Sersic 
index n as free parameter improves significantly the fitting 
and to find the best fitting structural parameters. The best 
fitting Sersic index n obtained for 6 out of the 9 galaxies 
is in the range 2.8 < n < 4.5. This result confirms that all 
of them are dominated by a bulge component. It is worth 
noting that the addition of the free parameter n does not 
significantly improve the fitting of these 6 galaxies with re- 
spect to the de Vaucouleurs profile, i.e. it is not statistically 
required as shown by the reduced we obtained. Thus, we 
can consider these galaxies morphologically as pure elliptical 
galaxies. For the remaining 3 galaxies, S7F5_254, S7F5_45 
and S2F1_443, a Sersic index n < 2.5 provides a better fit- 
ting to the data suggesting that a non negligible fraction, up 



to 30% of the light can be distributed on a disk component 
(see Saglia et al. 1997; van Dokkum et al. 1998). Given the 
spectrophotometric early-type nature of these galaxies (see 
Longhetti et al. 2005), we conclude that they probably are 
SO/Sa galaxies. 

Table 2 reports, for each galaxy, the total magnitude 
in the F160W filter and the morphological parameters we 
derived assuming the de Vaucouleurs and the Sersic pro- 
files. The reduced relevant to each best-fitting profile is 
also reported. The values of the effective radii (re and Re) 
are corrected for the small underestimate we found from 
our simulations (see next sub-section). The mean surface 
brightnesses {/i)f ^'''^^ and (p)^) have been derived from 
the corrected radii. 

Fig. 2 presents the surface brightness (SB) profile of the 
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Table 2. Morphological pajrameters of the galaxies. For each galaxy the first line reports the results obtained with the Sersic profile 
while the second line reports data obtained with n = 4 (de Vaucouleurs profile). Magnitudes are in the Vega system. F160W magnitudes 
can be transformed into AB mag by adding 1.36. The effective radii obtained with the Sersic and the de Vaucouleurs profiles have been 
corrected by about 0.07" and 0.025" respectively, a<;cording to the results of our simulations (see §3.2). The surface brightnesses have 
been derived from the corrected radii. 
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9 galaxies in the F160W band. The observed profile (filled 
points) extracted along the galaxy major axis is compared 
to the best-fitting de Vaucouleurs (solid/red line) and Sersic 
(dashed/green line) model profiles. The ellipse task within 
the IRAF package has been used to extract both the profiles 
on real and model images, adopting the best fitting param- 
eters (ellipticity and orientation) provided by Galfit. 

3.2 Simulations 

To assess the robustness of the results of the fitting to the 
profile of our galaxies, we applied the same fitting procedure 
to a set of simulated galaxies inserted in the real background. 
The main questions we want to address are i) what is the 
Sersic index we would obtain by fitting a Sersic profile to a 
pure n = 4 profile galaxy if placed at 2 ~ 1.4, sampled at the 
angular resolution of 0.075 arcsec (the pixel size of the NIC2- 
camera) and detected with a typical S/N> 20 and ii) which 
is the accuracy in recovering its effective radius. To this end, 
we generated with Galfit a set of 100 galaxies described by 
a de Vaucouleurs profile with axial ratio b/a and position an- 
gle PA randomly assigned in the ranges 0.4 < b/a < 1 and 



0<PA<180 dcg respectively. Magnitudes F160Wi„ and ef- 
fective radii rejn have been assigned randomly in the ranges 
19.5< FlGOWin <20.5 and 0.2 < r^,™ < 0.5 arcsec (corre- 
sponding to 1.7-4 kpc at 2: ~ 1.4) respectively. The simulated 
galaxies have been convolved with the HST-NIC2 PSF de- 
scribed in the previous section and embedded in the real 
background. The background image has been constructed 
by mosaicking different portions of real images devoid of 
sources. The sigma and mask images have been similarly 
derived from the real rms and mask images. We performed 
the Galfit fitting assuming the generalised Sersic profile 
(eq. 1) to the 100 galaxies. As first result, we find that they 
are best fitted by a Sersic profile with index 2.8 < Ufu < 3.9 
(< nfit >^ 3.3±0.2). It is worth noting that the value of the 
fitted Sersic index systematically decreases with fainter mag- 
nitudes, i.e. with lower S/N. This is shown in Fig. 3 where 
the value of Ufu is plotted versus the input magnitude of 
the simulated galaxies. This is possibly due to the fact that 
lower S/N values imply the missing of a higher fraction of 
the tail of a galaxy profile which would drop abruptly. This 
would favour lower values of n in the fitting. 

We then fit the simulated galaxies with the n = 4 pro- 
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Figure 2. In the upper panels of each figure, the surface brightness in the F160W band measured along the galaxy major axis of each 
galaxy (data points) is compared with the de Vaucouleurs (n = 4, solid/red line) and the Sersic [nfn, dashed/green line) profiles resulting 
from the fitting. The NIC2 pixel scale is 0.075 arcsec/pixel. In the lower panels the residuals of the fitting obtained as the difference 
between the data points and the models are shown. The formal la errorbar on the data points are shown for comparison. 



file and we studied the behaviour of the resulting r^ju with 
respect to the input re, in values for both the two profile fit- 
ting functions. In Figure 4 we plot the effective radius r^ju 
versus the input value re,m in the case of de Vaucouleurs 
and Sersic profiles. In the case of the n = 4 fitting profile 
the input effective radius of the galaxy is on average under- 
estimated by Arf ~ 0.025 arcsec (~ 10%; left panel 
of Fig. 4) while in the case of the Sersic fitting profile the 
mean underestimate is Arf'^''^'^ ~ 0.07" (~ 25%; Fig. 4 right 
panel). It is worth noting that this underestimate is smaller 
than the pixel size and apparently independent of input pa- 
rameters of the simulated galaxies (magnitude, axial ratio, 
size and PA). 



Finally, we exploited these simulations to assess a reli- 
able magnitude estimator for the galaxies in our images. To 
this end, we compared the SExtractor MAGJ3EST magni- 
tude and the Galfit magcai magnitude obtained for the 
100 simulated galaxies to the input FlGOWin magnitude. In 
Fig. 5 the input F160Wi„ magnitude is plotted versus the 
MAG_BEST magnitude (left panel) and versus the Galfit 
magnitude (right panel). While magcai is consistent with 
F160Wi„, MAG3EST is systematically fainter. The offset 
is about 0.2 mag with a rather small scatter of ±0.03 mag. 
On the basis of this result we used the Galfit magnitude as 
total magnitude of our real galaxies. 

The results obtained from the simulations confirm our 
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Figure 3. Best-fitting Sersic index Ufn obtained for the 100 sim- 
ulated n = 4 galaxies as a function of input magnitude. 
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Figure 4. Comparison between the effective radius of the sim- 
ulated galaxies (r^ i„) and the effective radius J^ejit obtained 
through the fitting with the de Vaucouleurs profile (left panel) 
and the Sersic profile (right panel). 



previous claims that the 6 galaxies with nju > 2.8 can be 
considered pure ellipticals. Only lower values of Ufu can be 
indicative of a deviation from a pure de Vaucouleurs profile 
suggesting the presence of a disk component in addition to a 
dominant bulge. Besides this, we also conclude that Galf it 
tends to underestimates the effective radius re in our im- 
ages by about 0.07 arcsec in the case of Sersic profile and 
by about 0.025" in the case of de Vaucouleurs profile. This 
underestimate has been taken into account in the derivation 



Figure 5. Upper panels - Comparison between the SExtractor 
(left) MAG_BEST magnitude and the Galf it magonj magnitude 
(right) with the input F160Wi„ magnitude of the 100 simulated 
galaxies. Lower panels - The difference A^ag = F^GOWi„ — mag 
is plotted as a function of MAG_BEST (left) and of magcal 
(right). 



of the mean surface brightnesses and will be considered in 
the following analysis. 



4 RESULTS 

The first result of the morphological analysis of our sample 
of early-type galaxies is that all of them are bulge dom- 
inated galaxies. In particular for 6 of them the best fit- 
ting Sersic profile does not improve the fit with respect to 
the de Vaucouleurs profile providing an index in the range 
2.8 < Ufit < 4.5. These galaxies can be considered pure 
elliptical galaxies. For the remaining 3 gala:xies (S7F5_254, 
S7F5.45 and S2F1_443) the values 1.9 < Ufu < 2.3 fit bet- 
ter the data suggesting the presence of a disk component 
possibly contributing up to 30% to the total light of the 
galaxy. 

A second piece of evidence is that the values of le 
obtained when the Sersic light profile is adopted to fit 
the galaxies light distribution arc systematically lower than 
those obtained with the de Vaucouleurs light profile. This ef- 
fect was expected also on the basis of the simulations show- 
ing that both the fitting profiles provide effective radii lower 
than the input ones. The underestimate is larger for Sersic 
profile fitting. We expect that this underestimate we observe 
at 2: ~ 1.4 would decrease with decreasing redshifts thanks 
to the higher spatial resolution and the larger radii at which 
galaxies can be sampled. This should be taken into account 
when the local morphological relations are compared with 
those at z > 1. 

Before doing this comparison, we recall the expected 
behaviour of the two fundamental morphological parame- 
ters used to describe the galaxies, i.e. Re [kpc] and (/Lt)e 
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[mag/arcscc'^]. As far as the first one, Re, it is usually as- 
sumed that it does not evolve and that it is a cliaractcristic 
scale of the early-type galaxies depending only from their 
mass. Prom this assumption, the expected evolution as func- 
tion of rcdshift of the mean surface brightness within Re can 
be deduced: 



{flx)e =Ma(2) + 5 log(Re) + 38.57 



(3) 



where M\{z) is the rest-frame absolute magnitude of the 
galaxy as function of z at fixed A. The expected surface 
brightness evolution of a galaxy, apart from the cosmological 
dimming, is due to the luminosity evolution expressed by its 
absolute magnitude at the fixed rest-frame wavelength. 

In the following section, we will compare the morpho- 
logical parameters measured on our sample of early-type 
galaxies at a ~ 1.4 with those of the local early-types. The 
comparison will be carried out on the basis of the Kormendy 
scaling relation involving Re and {fi}e- 



5 THE KORMENDY RELATION 

The Kormendy Relation (Kormendy 1977; KR hereafter) 
is a linear scaling relation between the logarithm of the 
effective radius Re [kpc] and the surface brightness /ie 
[mag/arcsec^] measured at Re: 



fj.e = a + /31og(Re) 



(4) 



The early-type galaxies both in field and in clusters follow 
this tight relation with a fixed slope /? ~ 3 out to a: ~ 1 (di 
Serego Ahghieri et al. 2005; Reda et al. 2004; La Barbera 
et al. 2003, 2004; Ziegler et al. 1999; Schade et al. 1996; 
Hamabe & Kormendy 1987). On the other hand, the zero 
point a is found to vary with the redshift of the galaxies 
and it reflects their luminous evolution over the time. Since 
the value of a strictly depends on the photometric band 
and system selected to derive the morphological parameters, 
its value needs to be transformed into that of a common 
rest-frame wavelength when comparisons at different z are 
performed. Furthermore, the KR is sometimes expressed in 
terms of the mean surface brightness (/x)e measured inside 
Re instead of /Xe, and the zero point of the two formalisms 
is different: 



(M)e = a'+/31og(Re). 



(5) 



If the projected light distribution for an elliptical galaxy is 
well described by the de Vaucouleurs profile, then 

q' = Q - 1.4 (6) 

and the more general relation between a and a' can be found 
in Graham & CoUess (1997). In the following we choose to 
adopt the KR expressed as in the formula (5), that is we 
consider the mean surface brightness (/x)e measured inside 
Re. Since our galaxies arc at 2 ~ 1.5, the observed F160W 
band (A = 16030A) corresponds to the rest-frame 6400A, 
that is about the rest-frame R-band. Thus, the comparison 
of our data with those of galaxies at lower redshift will be 
done on a KR derived in the standard Cousin R band. The 
different k-corrections applied to the various samples and 
the way we derived them are described in Appendix A. As 
far as the slope is concerned, we adopt (3 = 2.92 ± 0.08 as 
reported by the recent work by La Barbera et al. (2003). 



The zero point at z ~ can be derived from La Barbera 
et al. (2003) adopting the a value for the galaxies in the 
Coma cluster, a = 18.68. Morphological parameters of the 
galaxies in the Coma cluster come from j0rgensen, Pranx 
& Kjffirgaard (1995) and are measured in the Gunn r band. 
We transform a{r) into a{R) by adding 0.05 (due to slightly 
different filters) and -0.43 (due to the different photometric 
system of the gunn magnitudes): 



{n)e = 18.30 2.92 log(Re 



(7) 



If we consider zcoma = 0.024, for the cosmological dimming 
we derive a = 18.20 at « = 0. 

The zero point at 2: ~ can be also derived by the 
original Hamabe & Kormendy V-band relation for early- 
type galaxies in the local universe: 



He = 19.48-1- 2.92 log(Re) 



(8) 



Assuming a typical colour (V-R)=0.65 we derive the corre- 
sponding R-band a = 18.83 and by means of the equation 
(6): 



(M)e = 17.43-1- 2.92 log(Re). 



(9) 



It is worth noting that Hamabe and Kormendy possibly de- 
rive this relation assuming Ho — 50 km s~^ Mpc^^, even 
if not explicitly stated. If this was the case, their relation 
scaled to Ho — 70 km s~^ Mpc^^ would agree with the one 
derived by La Baxbera et al. (eq. 7). 

Fig. 6 report the KR expressed in the equation (7) (solid 
black line) and in eq. (9) (shaded green strip), this latter 
scaled to Hq — 70 km s^^ Mpc~^. The uncertainty of the 
KR by La Barbera et al. (2003) is shown as a strip limited 
by short-dashed black lines and the upper limit results co- 
incident with the original (uncorrected) KR by Hamabe & 
Kormendy (1987). The grey (red) dot-dashed line has been 
derived by the KR reported by La Barbera et al. (2003) 
for the cluster EIS0048 at z = 0.64. In the figure, our data 
are reported as laige filled circles. An open circle marks the 
three galaxies having n/a < 2.5 (S2F1_254, S2F1_443 and 
S7F5_45). In Fig. 6, the (/i)e and Re values of our galaxies 
have been measured through the Sersic profile fitting. The 
rest-frame surface brightness in the R-band, (/i)e, has been 
obtained from the observed F160W-band apparent magni- 
tudes. Since for each of our galaxies we have a detailed de- 
scription of its rest-frame optical spectral shape, we accu- 
rately derived the fc-correction: 

fc=mag(F160W;2; = Zgai)-ma,g{R,z = 0) 

where the two magnitudes are derived from the spectral en- 
ergy distribution reproducing that of each galaxy. ^ 

In Fig. 6 we also report data from literature resulting 
from both Sersic and dc Vaucouleurs profile fitting. Small 
filled symbols are the values relevant to the galaxies belong- 
ing to the samples of Reda et al. (2004, squares), Di Serego 
Alighieri et al. (2005, circles) and Moriondo et al. (2000, 

^ The uncertainty affecting this correction is mucli smaller tiian 
the photometric errors affecting the measure of (/ie) itself. This is 
due to the fact that the observed band at z = Zg^i is almost coin- 
cident with the rest-frame R band. Moreover, the sample of spec- 
tral templates used to calculate this fc-correction is constrained 
to few and similar spectra by the multi-band photometry and by 
the near-IR spectra available for all the galaxies (see §2). 
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Figure 6. Mean surface brightness (/i)e versus effective radius Re [kpc]. All the data are corrected for the cosmological dimming (l + z)'*. 
Our data (large filled circles) are compared with those of other authors. The large open circle superimposed on the filled one marks 
the three galaxies of our sample with rij-^t < 2.5 (S2F1_254, S2F1_443 and S7F5_45). The small filled symbols represent the data from 
Moriondo et al. (2000, MoriOO), Di Serego Alighieri et al. (2005, DiSe05) and Reda et al. (2004, Reda04). The open symbols represent 
the mean values of the surface brightness of field galaxies in the different redshift bins considered as derived from the data of Gebhardt 
et al. (2003, GebOS), Treu et al. (2005, Treu05) and Ferreras et al. (2005, Ferr05). The size of the symbols decreases with increasing 
redshift. The errorbar on the left marked with the open symbols represent the typical scatter in surface brightness within each bin. The 
solid line represents the KR at 2 ^ found by La Barbera et al. (2003) as expressed in eq. (7). The short-dashed lines represent the 
itlcr dispersion of the relation. The shaded (green) strip is the KR at z ^ of eq. (9) (scaled to Ho = 70 km s~^) found by Hamabe & 
Kormendy (1987). The long-dot-dashed (red) line is the KR at z ~ 0.64 by La Barbera et al. (2003). 



triangles). As 2: ~ reference, we report the sample of lo- 
cal galaxies of Reda et al. (2004) and derived assuming the 
de Vaucouleurs profile in the galaxy fitting. The morpho- 
logical parameters of the sample of Di Serego Alighieri et 
al. (2005) have been obtained by assuming a general Sersic 
profile. Their data have been transformed from B-band to 



R-band (see Appendix A). The sample of Moriondo et al. 
(2000) includes 6 galaxies at z ~ 1.3, 4 out of which have 
been observed with the same setup of our sample of galax- 
ies (HST-hNICMOS, F160W filter) while the remaining 2 
have been imaged in the HST-I-WFPC2 F814 fiher (see Ap- 
pendix A). As in the case of Di Serego Alighieri et al. (2005) 
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morphological parameters have been estimated by assuming 
general Sersic profiles. 

The open symbols represent the median values of the 
surface brightness of field galaxies in different redshift inter- 
vals as derived from the samples of Gcbhardt ct al. (2003, 
triangles), Treu et al. (2005; pentagon) and Fcrrcras et al. 
(2005; stars). The size of the symbols decreases with increas- 
ing redshift. The sample of Gebhardt et al. (2003) includes 
21 early-types with spectroscopic redshift and HST imag- 
ing in the F814W filter. They derive the cff'cctivc radii from 
a de Vaucouleurs fitting (Simard et al. 2002). The sample 
spans a very narrow range of Re, from 2.9 to 3.9 kpc. We 
divided this sample in the throe redshift bins 0.28 < z < 0.5 
(6 gal.) 0.5 < z < 0.8 (7 gal.), 0.8 < z < 1.0 (8 gal.) 
and we derived the surface brightness in the rest-frame R- 
band from the F814W apparent magnitude (see Appendix 
A). The sample of Treu et al. (2005) includes ~ 90 galaxies 
with morphological types (E) and 1 (E/SO), according to 
their classification. The morphological parameters have been 
obtained from a de Vaucouleurs profile fitting to HST- ACS 
images. We considered the four redshift bins 0.1 < z < 0.3, 
0.3 < z < 0.5, 0.5 < z < 0.7 and 0.7 < z < 1.0 and we 
derived {fJ-e)^' from the apparent magnitude in the F814W. 
The sample of Ferreras et al. (2005) includes ~ 70 galaxies 
with spectroscopic redshift in the intervals 0.2 < z < 0.5, 
Q.5< z < 0.7, 0.7 < 2 < 0.9 and 0.9 < z < 1.1. We selected 
them according to the photometric type tphot < 1.5. We 
derived the rest-frame R-band surface brightness from the 
F775W apparent magnitude. It is worth noting that the se- 
lection criteria used by Ferreras et al. to construct this sam- 
ple are quite different from those used for the other samples. 
Indeed, they classify the galaxies on the basis of the best- 
fitting template and then they exclude from the final sample 
those galaocies that are not consistent with an evolution into 
the Kormendy relation. 

Besides these data sets, a morphological study of early- 
type galaxies at 2: 1 has been made also by van Dokkum 
et al. (2003) and by van der Wei et al. (2004), represented 
in Fig. 5 by the starry symbols. The sample of van Dokkum 
et al. (2003) contains three galaxies selected from the clus- 
ter RDCS J0848+4453 at « = 1.27 observed with HST- 
NICMOS in the F160W filter. The sample of van der Wei 
ct al. (2004) is composed by 6 galaxies, 4 at z ~ 1 and 2 
at z = 0.7, observed with HST-ACS F850LP and F775W 
respectively. It should be noted that the effective radii they 
report in Tab. 1 are in arcscc and not in kpc as stated in 
the header of their table. In both samples the morphologi- 
cal parameters have been derived assuming de Vaucouleurs 
light profiles. 

5.1 Constraining the evolution of early- types 
galEixies 

Fig. 6 shows that the higher the redshift the more luminous 
are the galaxies at fixed Re. It is worth noting that all the 
data have been corrected for the cosmological dimming fac- 
tor (1 -f- z)"*. Moreover, the rest-frame wavelength sampled 
6500A) is weakly affected by episodes of recent star for- 
mation while it traces very well the bulk of the stellar mass of 
the galax;ies. Thus, the observed trend qualitatively reflects 
the luminosity evolution Mx{z) that galaxies arc expected 
to undergo with redshift. The position of the grey (red) dot- 



dashed line reveals that the increase of (/i)e and thus of the 
luminosity of early-types at z ~ 0.6 is at least 0.5-0.8 mag 
respect to 2 = for the same value Re. Approaching z ^ 1 
(see the samples of Di Serego Alighieri et al. 2005, Moriondo 
et al. 2000, Gebhardt et al. 2003 and Ferreras et al. 2005), 
a further increase in luminosity of 0.5-0.8 mag is evident. 
Our sample at 2: ~ 1.4 suggests a further increase of 0.5-0.8 
mag with respect to the sample at 2 ~ 1. By comparing our 
data with the KR of local ellipticals described by eq. (7) we 
derived the mean difference (A(/i)e) between the effective 
surface brightness of our early- type galaxies at z ~ 1.4 and 
the one of local early- types having the same Re. We find 
(A(/i)e) ~ 2.5 ± 0.5 mag, i.e. early-type galaxies at z ~ 1.4 
would occupy the KR at z ~ if their R-band luminosity 
decreases by about 2.5 mag at constant Re in the last 9 
Gyrs (from z ~ 1.4 to z ~ 0). It is also worth noting that, 
if we consider only the 6 pure elliptical galaxies of our sam- 
ple this estimate would slightly increase since the remaining 
galaxies show slightly smaller A(/i)e. This value is rather 
insensitive to the age of the galaxies of our sample as we 
have verified by considering separately the older from the 
younger galaxies. To verify if the evolution of (/i}e observed 
at < z < 1.5 is consistent with the luminosity evolution 
expected for early-type galaxies, we considered the maxi- 
mum luminosity dimming that a galaxy can experience in 
this redshift range. To this end we modeled a pure passive 
evolution with a Simple Stellar Population (SSP, Bruzual & 
Chariot 2003) 2.5 Gyr old at z ~ 1.5 (the mean age of our 
sample of early-types, see Tab. 1) and we derived the lumi- 
nosity evolution in the R-band at z ~ 0. We estimate that 
the maximum dimming that a galaxy with a stellar popula- 
tion 2.5 Gyr old can experience in 9 Gyr is about 1.6 mag. 
This value is very weakly sensitive to the assumed mean age. 
Indeed, it varies by ±0.1 mag for ages 2.5 ± 0.5 Gyr. Thus, 
our galaxies show an evolution of (/L()e which exceeds ~2.5 
times (~1 mag) the one expected in the case of luminosity 
evolution. This discrepancy cannot be simply recovered. A 
possibility is to hypothesize the evolution of the other quan- 
tity affecting (/i)e, i.e. the effective radius Re. If this is the 
case, the effective radius Re of early-types should increase 
at least by a factor ~ 1.5 from z ~ 1.5 to z ~ to account 
for this discrepancy. 

In the next subsection, we will follow a quantitative 
approach to the evolution of the KR and of the luminosity 
of the early-type galaxies as function of z. 

5.2 The evolution of the zero point of the 
Kormendy Relation 

It may be useful to recall that the zero point a of the KR 
represents at any redshift the surface brightness of galax- 
ies whose effective radius is Re = 1 kpc. The decrease (in 
magnitudes) of a at fixed A with 2 is a direct measure of the 
increase of the surface brightness of the galaxies with Re = 1 
kpc. Even in the case of a fixed slope of the KR with z, the 
increase of (/i)eiHe=i can be caused by two factors: 1) the 
increase of the absolute luminosity at fixed wavelength of 
galaxies with Re = 1 kpc assuming that Re does not evolve 
with z; 2) the decrease of Re with z of galaxies of fixed mass. 
It is generally assumed that only the first factor out of the 
two mentioned above is effectively active in the determina^ 
tion of the evolution of a in the KR. However, we will see 
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Table 3 . Zero point a of the KR derived from samples at z ~ 1 for 
fS = 2.92. The values in brackets have been derived considering 
only the 6 galaxies of our sample with Sersic index n > 2.8. 



Sample 



a(deVauc) 



a(Scr) 



OUR 

DiSeOS 

MoriOO 

DiSe05 + MoriOO 
Geb03 (z) ~ 0.85 
TreuOS (z) ~ 0.83 
FerrOS (z) ~ 1.0 



15.67(15.35)± 0.45 



16.41±0.4 
17.20±0.1 
16.52±0.1 



15.53{15.37)± 0.47 
16.64± 0.44 
16.09± 0.16 
16.50 ± 0.42 



that it could be necessary to hypothesize also the evolution 
of Re with z. 

In Fig. 7, we report the values of a as function of 
{1 + z) measured in the common rest-frame R band and 
derived from the samples described in the previous subsec- 
tion. For all the samples, we assume a fixed slope of the KR 
(/3 = 2.92). All the values have been already corrected for 
the cosmological dimming (1 -I- z)'^, and thus they directly 
describe the evolution of the zero point of the KR with z due 
to the luminosity and/or scale evolution of the galaxies. The 
black filled squares at z ^ 0.64 come from La Barbera et al. 
(2003), the small filled circle at z = 1 has been derived from 
the data of Di Serego Alighieri et al. (2005), the filled tri- 
angle from the sample of Moriondo et al. (2000). The open 
symbols come from the samples of Gebhardt et al. (2003, 
triangles), Treu et al. (2005, pentagon) and Ferreras et al. 
(2005, stars). The largest filled circles reported at z = 1.4 
represent our results in the case of Sersic profile fitting (open 
circle) and de Vaucouleurs profile (filled circle). In Table 3 
the values of a we derived at z ~ 1 are reported. 

In Fig. 7, the expected evolution of a due to the lumi- 
nosity evolution of galaxies with fixed mass (i.e., due to the 
decrease of M/L) for some models is superimposed on the 
data. In the upper panel all the model lines are derived as- 
suming a fixed star formation history (SFH, r — 0.6 Gyr) at 
solar metallicity, and they differ for the formation redshift 
Zf = 2, 3, 4, 6 from the top to the bottom. In particular, the 
black solid line represents the model with Zf — 4. In the 
bottom panel, all the models refer to Zf — 4 and they dif- 
fer for the SF history assumed, where r = 2.0,0.6,0.1 Gyr 
from the top to the bottom. The black solid line represents 
the same model as in the upper panel. All the models are 
obtained with the BC03 code, and they assume Chabrier 
(2003) IMF. We verified that the choice of diflferent IMFs 
does not change the results. 

Fig. 7 shows that data at z ^ 1 are only marginally 
compatible with the expected evolution of the luminosity 
of the galaxies. Indeed, they almost match only lines corre- 
sponding to models which include star formation activity at 
1 < 2 < 2, i.e. 2/^2 for r = 0.6 Gyr (upper panel) or r = 2 
Gyr for Zf ^ 4 (lower panel). On the other hand, the spec- 
trophotometric analysis of the whole SED of at least some 
of these samples of galaxies, from UV to near IR, shows 
that the bulk of their stellar content must be older than 
some Gyr. Furthermore, many of the samples represented in 
Fig. 7 have been selected on the basis of red optical-near IR 
colours (e.g., R-K>5) which are not compatible with stellar 
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Figure 7. Zero points a of the Kormendy relations derived by the 
various samples as a function of their redshift (/3 = 2.92). Large 
filled (empty) circle represent the zero point we obtain from our 
sample in the case of de Vaucouleurs (Sersic) profile. The lines 
show the expected evolution of a due to the luminosity evolution 
of galaxies, i.e. due to the evolution of their M/L ratio for different 
models. Upper panel: models refer to the same SFH r = 0.6 Gyr 
starting at the redshifts of formation (from the top to the bottom) 
Zf = 2,3,4,6. Lower panel: models assume the same Zf = A but 
refer to different SFHs, from the top, r = 2.0, 0.6, 0.1 Gyr. 



contents younger than 1 Gyr on average. Thus, if we try 
to interpret Fig. 7 on the basis of the luminosity evolution 
of galaxies we hardly reconcile models and data. A possible 
way to solve this problem is to assume that galaxies evolve 
also in their morphological scale, that is to admit that the 
effective radius of a galaxy of fixed mass increases from the 
epoch of its formation toward 2; = 0. As we have already 
recalled, the zero point a of the KR represents the surface 
brightness of the galaxies with Re = 1 kpc, and from eq. 2 it 
is expected to evolve as MA(z)[ji^=i] once the cosmological 
dimming has been taken into account: 



(MA)e|fle = l ~ MA(z)[fl^ = l] 

It is usually assumed that Ma(2)[h^=i] varies only as a con- 
sequence of the variation of the stellar mass to light ratio: 

MA(2)[fl^=i] « -2.51og[(LA/M)(2)]-2.51og(Mass)[«^=i](ll) 

where only Lx/M is function of 2, while the stellar mass 
content of the galaxies with Re = 1 kpc is constant. 
Assuming evolution of Re from 2/ to 2 = implies that, at 
each redshift, galaxies of a given Re have variable mass. The 
expected change in the value of a must contain an additive 
increase with 2 due to the increase of mass of the galaxies 
with fixed scale, that is Mass[/j^^i] becomes Mass[^^=i] (2) 
function of z. Fig. 8 shows an example of the expected evo- 
lution of a with 2 in the assumption that the evolution of Re 
stops for 2 < 0.5. Points and lines represent the same data 
and models as in Fig. 7, but for 2 > 0.5 the value of the zero 
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Figure 8. Same as Fig. 8 but models include the evolution of the 
effective radius Re, as described in §5.2. 

point of the KR increases by A(a) = —2.5 log(0.5 + z), that 
corresponds to assume: 

Re(^>0.5) = i^^. (12) 

for a fixed stellar mass of a galaxy. In this case, the mod- 
els which describe the spectrophotometric properties of the 
galaxies, i.e. ~ 4 and SF histories with short time scales 
(t < IGyr) can reproduce also the observed behaviour of 
a. The evolution of the size of the spheroid-like objects has 
been recently hypothesized also by other authors, at least 
for the most massive galaxies. For instance, Trujillo et al. 
(2006b) find strong hints that the most massive early-type 
galaxies were ~ 4 times smaller at z > 1.4 than at z = 0. If 
we consider z — 1.4— 1.5 in eq. (12) we obtain a lower (even 
if comparable) evolution. It is important to point out that 
the increase of the stellar mass at 2 > 0.5, needed to explain 
the observed evolution of the zero point of the KR, can be 
hardly ascribed to on-going SF activity at 0.5 < z < 1.5. 
Indeed, from one hand the observed SED of the early-type 
galaxies at 2~ 1.0 — 1.5 is not compatible with recent or 
on-going SF (Longhetti et al. 2005) and, on the other hand, 
our sample at 2 ~ 1.4 is composed by galaxies with stellar 
masses Metar ^ 3 X 10^^ M© which cannot form more stars 
at z < 1.4 to be compatible with the local mass function 
of the early-type galaxies. Thus, the possible solution is to 
hypothesise that Re of early-type galaxies increases from at 
least z ~ 1.5 toward z = and that this behaviour is not 
related to a growth of their stellar mass by star formation. 
A possible evolution of the effective radius is suggested also 
from the observation that the values of Re measured for the 
galaxies of our sample are smaller on average than the values 
measured on galaxies at lower redshift (see Fig. 6) in spite 
of the high-mass of our galaxies. Finally, it is worth noting 
that the assumed scale evolution of galaxies corresponding 
to the increase of Re from their formation toward z = does 



not imply any change in the slope of the KR. It can be seen 
as single galaxies during their ageing move along the KR, 
causing a further change in its zero point but not a change 
in its slope. 



6 SUMMARY AND CONCLUSIONS 

We presented the morphological analysis of a sample of 9 
massive (Mstar > 10^^ M0) galaxies spectroscopically clas- 
sified as early-types at 1.2 < z < 1.7. The analysis is based 
on HST imaging carried out with the NICMOS-NIC2 cam- 
era (0.075 arcsec/pixel) in the F160W filter (A ~ 16000 A), 
sampling the rest-frame R-band. We studied their profile 
and we derived their effective radius re and their mean sur- 
face brightness (/i)e within Re by means of the Galf it soft- 
ware (v. 2.0.3; Peng et al. 2002). From the surface bright- 
ness fitting analysis performed we find that the 9 early-type 
galaxies analysed are bulge dominated in agreement with 
the spectral classification. Six of them are very well fitted 
by the pure de Vaucouleurs profile (n = 4) and the addi- 
tion of the index n as a free parameter (providing values 
2.8 < Ufit < 4.5) is not statistically required by the fitting. 
Thus, these 6 galaxies can be considered pure ellipticals. For 
the remaining 3 galaxies the values 1.9 < rifu < 2.3 suggest 
the presence of a disk component possibly contributing up 
to 30% to the total light of the galaxy suggesting that they 
are SO galaxies. The median effective radius of our sample 
of early-types is ~ 2.95 kpc and ~ 2.8 kpc in the case of de 
Vaucouleurs and Sersic profile respectively. 

By applying the same analysis to a set of simulated 
galaxies with intrinsically pure n — 4 profile embedded in 
the real background we find that they are fitted by values 
2.8 < Ufit < 3.9 confirming the previous morphological clas- 
sification. 

The data we have analysed allow us to exclude that the 
sizes measured are affected by morphological k-corrections 
as (possibly) in the case of blue and UV observations, as 
suggested by Daddi et al. (2005). Indeed, our observations 
sample the rest-frame R-band mainly sensitive to the old 
stellar population. Moreover, the high resolution (FWIIM~ 
0.1 arcsec) of our IIST-NIC2 images, corresponding to a 
spatial sampling of ~ 0.8 kpc at the redshift of our galaxies, 
allow us to exclude that the estimate of the effective radii is 
limited by the resolution. 

The analysis of the KR reported to a common rest- 
frame R-band of our data compared with other samples at 
lower redshift shows that the higher the redshift the more lu- 
minous are the galaxies at fixed Re, as qualitatively expected 
from the luminosity evolution of galaxies. However, we find 
that the mean surface brightness (/i)e of our sample should 
evolve by about 2.5 mag to z ~ to occupy the local KR im- 
plying a luminosity evolution of an order of magnitude in the 
last 9 Gyr. This luminosity evolution exceeds by a factor 2.5 
(1 mag) the pure passive evolution, the maximum evolution 
expected for an elliptical galaxy. This suggests that other 
parameters, possibly structural, may undergo an evolution. 

The observed effect mentioned above has been quanti- 
fied also by means of the evolution of the zero point a of the 
KR as function of z. The resulting evolution of a cannot be 
simply explained in terms of the evolution of the luminosity 
of the galaxies, at least at z > 0.8. Indeed, the strong de- 



The Kormendy relation at z ^1.5 13 



crease of a{z) at 2 > 0.8 — 1.0 requires to assume that the 
stellar content of the early-type galaxies at that redshift is 
too young (< 1 Gyr) and that they are still forming stars. 
But these assumptions axe not compatible with their global 
spectrophotometric properties which date their stellar con- 
tent as old (3> 1 Gyr). 

A way to reconcile models and observations is to hy- 
pothesise that mori)hological scale plays a role in the evo- 
lution of the surface brightness. We proposed a scheme in 
which the effective radius of a galaxy of fixed mass increases 
from the epoch of its formation toward « = 0. In partic- 
ular, if Re(2;) = ''^('J^osf ^ ' same models describing the 
spectrophotometric properties of the early-types observed at 
z > 1 are able to match the observed evolution of a with 
z. On the other hand, Waddington et al. (2002) present evi- 
dence that the smaller size of their 2 radio galaxies at z = 1.5 
compared with the sample of more powerful radio galaxies 
at 2 ~ 0.8 by McLure & Duulop (2000) can be explained 
simply assuming a lower mass without any need of scale 
evolution. It is thus possible that the evolution of Re with 
z that we have proposed involves only galaxies with stellar 
meisses Mstar > 10^^ M©. The evolution in the galaxy sizes 
seems to be indeed dependent on the mass of the objects as 
pointed out by many authors (e.g. Treu et al. 2005; Trujillo 
et al. 2006a, 2006b). However, the zero points of the KR de- 
rived from the samples of Di Serego Alighieri et al. (2005), 
Moriondo et al. (2000), Gebhardt et al. (2003) and Ferreras 
et al. (2005) which contain less massive galaxies, present the 
same inconsistencies with models as those derived from our 
sample, even at a lesser extent. 

Obviously, many problems arises from the assumption 
of size evolution, and questions on the physical framework in 
which such an evolution can take place. All our galaxies have 
high stellar masses (M^tar- > lO^'^ M©) and they account for 
almost all the number density of local early-typo galaxies 
with comparable mass (Saracco et al. 2005). Thus, not all of 
them can grow further their mass from z ~ 1.5 to z = to be 
compatible with the local mass function of early-types, set- 
ting aside the process responsible of the growth. This implies 
that mergers (dry or wet) at 2 < 1.5 cannot play a major role 
in the evolution of high-mass galaxies in this redshift range, 
at variance with the findings of some authors (e.g. Bell et 
al. 2006). In fact, dry mergers possibly combined with AGN 
feedback were proposed to reconcile the observed properties 
of early-type galaxies, including the possible size evolution, 
with the hierarchical paradigm (Khochfar and Burkert 2003; 
Khochfar and Silk 2006; De Lucia et al. 2006; Hopkins et al. 
2006). Dry mergers should occur among gas-poor progen- 
itors of evolved stars without additional star formation to 
form massive early-type galaxies. Since dry mergers are less 
affected by dissipation mechanism than mergers involving 
large amounts of gas (wet mergers), the stellar components 
of the final products are less centrally contracted and, con- 
sequently, the optical sizes of dry mergers products should 
be larger than the wet mergers products. It is difficult to 
arrange this scenario to our results. Indeed, since not all of 
our galaxies can grow further their stellar mass, we can hy- 
pothesize that only a fraction (~30%) of them experience 
dry merger at 2 < 1.5 accounting for the local density of 
Mstar ~ 10^^ Mq galaxies. However, the remaining frac- 
tion of galaxies would still remain without a counterpart at 
2 ~ 0. 



If our galaxies have already completed their growth at 
2 ~ 1.5, dry mergers could be the mechanism responsible of 
their assembly at 2 > 1.5 provided that something happens 
at 2 < 1.5 to enlarge their apparent sizes. On the contrary, 
the hypothesis that our early-types arc the products of wet 
mergers, which would fit better their small sizes observed 
at 2 ~ 1.5, hardly fit their other properties. Indeed, the old 
age characterizing their stellar populations and the relevant 
formation redshift 2 > 4 (Longhetti et al. 2005) cannot be 
accounted for by wet mergers which would induce star for- 
mation episodes. 

Thus, we leave with the evidence in hand that high-mass 
early-type galaxies at 2 ~ 1.5 arc more compact then their 
local counterparts. If this is due to a structural evolution 
we should consider that smaller sizes imply higher densities 
and, consequently, higher velocity dispersions. In particular, 
since our galaxies are at least 1.5-2 times smaller than those 
at 2 ~ they should have velocity dispersions at least 1.2- 
1.5 times larger, i.e. larger than ~ 300 km/s. However, such 
high-mass and high-density systems are not seen in the lo- 
cal universe, at least with this high number density. Thus, 
in the hypothesis of a structural evolution something must 
happens at 2 < 1.5 to relax these systems and to match the 
observed size and velocity dispersion of the local population 
of high-mass early-types. On the basis of these considera- 
tions we believe that a key observational test would be the 
measure of the velocity dispersion of such galaxies and a 
spatial map of their kinematics. These data would allow us 
to unambiguously address the question whether or not the 
smaller sizes observed in our massive early-type galaxies are 
due to higher densities, i.e. to structural and dynamical evo- 
lution occurring after the completion of their assembly. 



APPENDIX A: K-CORRECTIONS APPLIED TO 
LOWER REDSHIFT SAMPLES 

All the k-corrections calculated throughout this paper are 
based on the latest version of the Bruzual & Chariot models 
(2003; BC03), and they are obtained averaging the values 
which can be derived assuming the Salpeter (1955) and the 
Chabrier (2003) Initial Mass Function (IMF). 

The original KR reported by La Barbera et al. (2003) 
for the cluster EIS0048 at 2 = 0.64 is in the I band, and 
the corresponding zero point has been transformed into the 
common rest-frame R-band by applying a k-correction of - 
0.02 mag assuming a typical spectral shape of a template 
at 2 = 0.6 with solar metallicity, 10 Gyr old, and built by 
means of an exponentially decaying SFR with a time scale 
r =0.6 Gyr. 

No corrections have been applied to the Reda et al. 
(2004) data since they provide the measured quantities in 
the R-band. 

Di Serego Alighieri et al. (2005) provide measure of re 
and (/i)e in the B band. We transformed their morphological 
parameters in our common rest-frame R-band by applying 
an average k-correction of -1.4, that corresponds to 2 ~ 1 
templates at solar metallicity, 4 Gyr old and obtained as- 
suming r = 0.6 Gyr. 

The sample of 6 galaxies at 2 ~ 1.3 by Moriondo, 
Cimatti & Daddi (2000) is composed of 4 galaxies ob- 
served with the same setup of our sample of galaxies 
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(HST+NICMOS, F160W filter) and 2 galaxies in the 
HST+WFPC2 F814 filter. The values of ^lABmag have been 
transformed into our standard common R-band rest-frame 
by applying k-corrections of +0.7 and -1.7 for the F160W 
filter and F814 filter respectively [z ~ 1.3 template, at solar 
metallicity, 4 Gyr old, t = 0.6 Gyr). Another further cor- 
rection applied to this set of data is due to the fact that the 
authors compute instead of (/i)e. We adopt the Grahams 
& CoUess (1997) relation to obtain this correction factor. 

Gebhardt et al. (2003) provide measure of Re, of the 
absolute magnitude M_g and of the {^)e in the rest-frame B 
band in their Tab. 1. We divided the sample in three redshift 
intervals with median redshift 0.4, 0.63 and 0.85. For each 
galaxy we derived (^)e in the R-band rest-frame using the 
apparent magnitude F814W and applying a k-correction of 
0.35 at 2 ~ 0.4, 0.15 at z ~ 0.63 and -0.25 at z ~ 0.85. We 
considered templates at solar metallicity with r = 0.6 Gyr 
with age 8, 7 and 5 Gyr respectively. 

Treu et al. (2005) provide a sample of more than 220 
galaxies with morphological type at 0.2 < z < 1.2. We se- 
lected those having morphological type (E) and 1 (E/SO) 
and divide them into four redshift intervals centered at the 
median redshift 0.25, 0.45, 0.6 and 0.85. We derived in 
the rest-frame R-band from the photometry in the F814W 
band. Also in this case, we applied to each galaxy a k- 
correction estimated at the median redshift of the relevant 
interval, namely 0.5, 0.3, 0.17 and -0.25 mag derived by tem- 
plates 10, 8, 7 and 5 Gyr old respectively (solar metallicity 
and r = 0.6 Gyr). 

The sample of Ferreras et al. (2005) has been distributed 
in the four redshift interval centered at 0.38, 0.62, 0.74 and 
1.0. We derived {/i)e in the rest-frame R-band from the pho- 
tometry in the F775W band. The k-corrections relevant to 
the redshift intervals are 0.35, 0.15, -0.04 -0.6 mag. 

The sample of van Dokkum et al. (2003) contains three 
galaxies selected from the cluster RDCS J0848-I-4453 at z = 
1.27 observed with HST-NICMOS in the F160W filter. The 
k-corrections to transform their data into the standard rest- 
frame R-band are provided by the authors. 

The sample of van der Wei et al. (2004) is composed by 
6 galaxies, 4 at z ~ 1 and 2 at z = 0.7, observed with HST- 
ACS F850LP and F775W respectively. Small corrections are 
needed to be applied to the photometric data of this sample 
to be transformed in the standard rest-frame R-band (i.e., 
4-0.05, corresponding to the average value derived from tem- 
plates at the 2 different redshift, at solar metallicity, with 
age between 4 and 5.5 Gyr, obtained assuming r = 0.6 Gyr). 
For both the samples, the surface brightness provided by the 
authors is instead of (^i)e. Since for both the samples the 
morphological parameters have been derived assuming de 
Vaucouleurs light profiles, we transformed into by 
means of eq. (6). 
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